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Abstract. A very classical subject in Commutative Algebra is the Invariant 
Theory of finite groups. In our work on 3-dimensional topology 1121 . we found 
certain examples of group actions on polynomial rings. When we tried to 
compute the invariant ring using Singular [6] or MAGMA pQ, it turned out 
that the existing algorithms did not suffice. 

We present here a new algorithm for the computation of secondary invari- 
ants, if primary invariants are given. Our benchmarks show that the imple- 
mentation of our algorithm in the library f invar of Singular [6] marks a 
dramatic improvement in the manageable problem size. A particular benefit 
of our algorithm is that the computation of irreducible secondary invariants 
does not involve the explicit computation of reducible secondary invariants, 
which may save resources. 

The implementation of our algorithm in Singular is for the non-modular 
case; however, the key theorem of our algorithm holds in the modular case as 
well and might be useful also there. 
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1. Introduction 

Let G be a finite group, linearly acting on a polynomial ring R with n variables 
over some field K. We denote the action of g € G on r € R by g.r G R. 

Let R G = {r e R: g.r = r, Vg £ G} be the invariant ring. Obviously, it is a 
sub- algebra of R, and one would like to compute generators for R G . We study here 
the non- modular case, i.e., the characteristic of K does not divide the order of 
G. Note that according to [5], algorithms for the non-modular case are useful also 
in the modular case. 

For any subset S C R, we denote by ((S)) C R the sub-algebra generated by 
S, and by (5) C R the ideal generated by S. It is well known [2] that there 
are n (the number of variables) algebraically independent homogeneous invariant 
polynomials P = {pi, ...,p n } C R such that R G is a finitely generated {{P))~ 
module. The elements of P are called primary invariants. Of course, they 
are not uniquely determined. There are various algorithms to compute primary 
invariants [8]. Since the primary invariants are algebraically independent, the sub- 
algebra ((P)} is isomorphic to a polynomial ring with n variables. It is called 
(homogeneous) Noetherian normalization of R G . 

Let S C R G be a minimal set of homogeneous ((P))-module generators of R G . 
The elements of S are called secondary invariants. Note that the number of 
secondary invariants depends on the degrees of the primary invariants. Hence, it is 
advisable to minimize the degrees of the primary invariants. Irreducible secondary 
invariants are those non-constant secondary invariants that can not be written as 
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a polynomial expression in the primary invariants and the other secondary invari- 
ants. The set of secondary invariants is not unique, even if one fixes the primary 
invariants. It is easy to see that one can choose secondary invariants so that all of 
them are power products of irreducible secondary invariants. 

The aim of this paper is to present a new algorithm for the computation of 
(irreducible) homogeneous secondary invariants, if homogeneous primary invariants 
P are given. The key theorem for our algorithm concerns Grobner bases and holds 
in arbitrary characteristic; however, the algorithm assumes that we are in the non- 
modular case. For simplicity, we even assume that K is of characteristic 0, but this 
is not crucial. 

The rest of this paper is organised as follows. In Section [2j we briefly expose 
our motivating examples arising in low-dimensional topology. In the Section [3] we 
recall the basic scheme for computing secondary invariants. In Section [4] we state 
our key result and formulate our new algorithm for the computation of (irreducible) 
secondary invariants. In Section [5] we provide some examples (partially inspired 
by our study of problems in low-dimensional topology) and compare the implemen- 
tation of our algorithm in Singular [6] with previously implemented algorithms 
in Singular by A. Heydtmann [7] respectively in Magma p] by A. Steel [TO] . 



2. Motivating examples 

The starting point of our work was the study of generalisations of Turaev- Viro in- 
variants [llj , [12] . These are homeomorphism invariants of compact 3-dimensional 
manifolds. Their construction is (with some simplifications) as follows. Let T be 
some finite set, and let T be a triangulation of a compact 3-manifold M. An 
^-"-colouring of T assigns to any edge of T an element of T . Tetrahedra have six 
edges. So, for any tetrahedron of T, an ^"-colouring of T gives rise to a six-tuple 
of colours, that is called 6j -symbol and denoted by | J \ c d | , for a, b, c, d,e, f E T . 
The equivalence classes of 6j-symbols with respect to tetrahedral symmetry are 
variables of some polynomial ring, R. The ring also contains one variable Wf for 
any / e T, called the weight of /. For any JF-colouring, we form the product 
over the weights of the coloured edges and over the 6j-symbols of the coloured 
tetrahedra of T. By summation over all possible ^-colourings of T, we obtain a 
polynomial TV(T) called the state sum of T . Due to the tetrahedral symmetry of 
the 6j-symbols, the state sum is well-defined. However, it depends on the choice 
of T rather than on the homeomorphism type of M . It was shown by V. Turaev 
and O. Viro [16] that an appropriate evaluation of the state sum (yield by the rep- 
resentation theory of Quantum Groups) is independent of the choice of T . This is 
called a Turaev- Viro invariant. 

In [IT] and [12], we define an ideal I C R, the Turaev-Viro ideal. We show 
that the coset tv(M) = TV(T) + 1 is independent of T, hence, a homeomorphism 
invariant of M. This generalises the classical Turaev-Viro invariants. By extensive 
computations, we show in [12j that these so-called ideal Turaev-Viro invariants are 
much stronger than the classical Turaev-Viro invariants. For this, it was necessary 
to compute Grobner bases of Turaev-Viro ideals. It turns out that different al- 
gorithms for the computation of Grobner bases differ widely in their performance. 
The algorithm slimgb in Singular [6] of M. Brickenstein [2] performs particularly 
well. 
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We obtain a lower bound for the number of tetrahedra of any triangulation of 
M, in terms of the minimal degree of polynomials in the coset tv{M). However, in 
our computations, the bound appears to be trivial [12] . There was some hope to 
improve the lower bound as follows, using computations of invariant rings. Let G be 
the symmetric group of T . In the obvious way, G acts on the tetrahedral symmetry 
classes of 6j-symbols and on the weights, and hence, on R. The G-action permutes 
the summands of the state sum. So, the state sum belongs to R G . Let I G = IC\R G . 
Obviously, if T is a triangulation of some compact 3-manifold M, then the coset 
tv G (M) = TV(T) + I G C R G is a homeomorphism invariant of M, and as such in 
fact equivalent to tv{M). However, since tv (M) C tv(M), there is some hope that 
the minimal degree of polynomials in the coset tv (M) is higher than in tv(M), 
which would provide stronger bounds for the number of tetrahedra. 

This is how we became interested in the computation of invariant rings. The 
existing implementations in Magma and Singular could not compute the sec- 
ondary invariants in several of our examples. This motivated us to develop a new 
algorithm for the computation of secondary invariants. It has been part of the 
f invar library of Singular [6] since release 3-0-2 (July 2006). Unfortunately, in 
our topological applications, we did not find an improvement of the lower bound 
for the number of tetrahedra. However, our new algorithm for the computation of 
secondary invariants certainly is of independent interest. 

3. Generalities on the computation of secondary invariants 

In the non- modular case, we can use the Reynolds operator Rey: R — > R G , 
which is defined by 

Rey(r) = ^ g.r 
1 1 sec 

for r G R. By construction, the restriction of the Reynolds operator to R is the 
identity. Let Bd C R G be the images under the Reynolds operator of all monomials 
of R of degree d. It is well known that one can find a system of homogeneous 
secondary invariants of degree d in Bd [15] . But how can one determine what 
elements of Bd are eligible as secondary invariants? 

Let So, S\, 52, Sd-i C R G be the homogeneous secondary invariants of de- 
gree 0, 1, 2, d — 1, respectively (we can take Sq — {1}), and let I Si c Si be the 
irreducible ones, for i = 1, d — 1. Let si, s m 6 R G be some homogeneous sec- 
ondary invariants of degree d. Let b € Bd- We can choose 6 as a new homogeneous 
secondary invariant of R G , if b is not contained in the ((P))-module generated by 
So U Si U • • • U Sd-i U {si, s m }. It is not difficult to show that this is the case if 
and only if b is not contained in the ideal (P U {si, s m }) C R; see [15] . 

Ideal membership can be tested using Grobner bases. For p 6 R and a finite 
subset Q C R, we denote the remainder of p by reduction modulo Q by rem(6; Q). 
The remainder is iteratively defined, depends on the choice of a monomial order, 
and in general depends on the order of the elements of Q. For a definition of 
remainder, of Grobner bases, and for a proof of the following classical result, we 
refer to [5] or [14] . 



Theorem 1. Let Q be a Grobner basis of (G) C R, and let p E R. Then, rem(p; Q) 
does not depend on the order of polynomials in Q , and we have rem(p; Q) — if and 
onlyifpe(G). □ 
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We thus obtain the following very basic algorithm for finding homogeneous sec- 
ondary invariants Sd of degree d, provided those of smaller degrees have been com- 
puted before. 
Basic Algorithm 

(1) Let Sd = 0- Let Q be a Grobner basis of (P). 

(2) For all b e B d : 

If b £ (P U Sd) (which is tested by reduction modulo Q) then replace Sd by 
S d U {&}; compute a Grobner basis of (P Li S d ) and replace Q with it. 

(3) Return Sd- 

There are several ways to improve this algorithm. One way is an application 
of Molien's Theorem [15], [8], [7]. We will not go into details here. Molien's 
Theorem allows to compute the number of secondary invariants of degree d. 
In other words, if in the above algorithm we got md secondary invariants, we can 
immediately break the loop in Step (2). 

We also would like to see which of the secondary invariants in Sd are irreducible, 
since these, together with P, generate R G as a sub-algebra of R. For that purpose, 
one forms all power products of degree d of elements of I Si U IS2 U • • • U ISd-i and 
chooses from them as many secondary invariants as possible (compare [7] or [5]). If 
there are further secondary invariants (which we know from computation of md), 
then one proceeds as above with Bd, and obtains all irreducible secondary invariants 
ISd of degree d. So, the algorithm is as follows. 
Refined Algorithm 

(1) Compute md- Let Sd — IS d — and let Q be a Grobner basis of (P). 

(2) For all power products b of degree d of elements of IS\ Li IS2 U • • • U ISd-i- 

(a) If b (PL) S d ) (which is tested using Q) then replace Sd by Sd, U {&}; 
compute a Grobner basis of (P Li Sd) and replace Q with it. 

(b) If \Sd\ — md then break and return (Sd,ISd)- 

(3) For all b 6 B d : 

(a) If b (P Li Sd) (which is tested using Q) then replace Sd by Sd U {&}, 
and ISd by ISdL){b}; compute a Grobner basis of [P Li Sd) and replace 
Q with it. 

(b) If ISj I = rrid then break and return (Sd,ISd)- 

Eventually, Sd contains homogeneous secondary invariants of degree d, and ISd 
contains the irreducible ones. In this form, the algorithm has been implemented 
in 1998 by A. Heydtmann [7] as the procedure secondary_charO of the library 
f invar of Singular. In Step (2), the ideal membership is tested by computing 
the remainder modulo some Grobner basis of the ideal. This ideal changes once a 
new secondary invariant has been found. So, the algorithm involves many Grobner 
basis computations. This is its main disadvantage and limits the applicability of 
the Basic and the Refined Algorithm. 

An alternative algorithm was proposed by Kemper and Steel (see jS], [TU] or [5]) 
and implemented in Magma pQ. Here, new secondary invariants are detected not 
by a general solution of the ideal membership problem but by Linear Algebra. 
This algorithm only involves one Grobner basis computation, namely for the ideal 
(P). But for computing some of the invariant rings that arise in our study of 
homeomorphism invariants of 3-dimensional manifolds [12] . this does not suffice 
either. 
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4. The New Algorithm 



The main feature of our new algorithm is that, after computing some (homoge- 
neous) Grobner basis of (P), we can directly write down a homogeneous Grobner 
basis up to degree d of {PUSd}, once a new secondary invariant of degree d has been 
found. We can do so whithout any lengthy computations (in contrast to [7]), and 
we also avoid to deal with huge systems of linear equations (in contrast to [10] , [8] , 
[3j). This allows to solve the ideal membership problem in a very quick way. We 
recall the notion of "homogeneous Grobner bases up to degree d" in the following 
paragraphs. At the end of the section, we provide our key theorem and formulate 
our new algorithm. 

For p £ R, let lm(p) the leading monomial of p, let lc(p) be the coefficient 
of lm(p) in p, and let lt(p) — lc{p)lm(p) be the leading term of p. The least 
common multiple is denoted by LCM(-, •). Now we can recall the definition of the 
S polynomial of p, q £ R: 



Obviously, the S'-polynomial of p and q belongs to the ideal (p, q) C R. The 
leading terms of p and q are canceling one another, so, the leading monomial of 
S(p, q) corresponds to monomials of p or q that are not leading. The following 
result can be found, e.g., in [5] or [T3]. 

Theorem 2 (Buchberger's Criterion). A set gi,...,g k £ R of polynomials is a 
Grobner basis of the ideal {g\, ...,gfc) C R if and only ifrem(S(gi,gj);gi, ...,3fe) = 
for all i,j = 1, ...,k. □ 

Buchberger's Criterion directly leads to Buchberger's algorithm for the construc- 
tion of a Grobner basis of an ideal: One starts with any generating set of the ideal. 
If the remainder modulo the generators of the S-polynomial of some pair of gen- 
erators does not vanish, then the remainder is added as a new generator. This 
will be repeated until all S'-polynomials reduce to 0; it can be shown that this will 
eventually be the case, after finitely many steps. 

Here, we are in a special situation: We work with homogeneous polynomials. It 
is easy to see that if p and q are homogeneous then so is S(p,q), and its degree 
is higher than the maximum of the degrees of p and q, unless lm(p) — lm(q). 
If p, g\,g%, ...,<7fe £ R are homogeneous then so is rem(p; gi, g^j. Moreover, 
either rem(p; g x , g k ) = or deg (rem(p; gi, g k )) = deg(p). For computing 
rem(p; gi, gfc), only those gi play a role with deg(<?i) < deg(p), for i = 1, 
It follows: If an ideal I C R is homogeneous (i.e., it can be generated by homoge- 
neous polynomials) then it has a Grobner basis of homogeneous polynomials. Such 
a Grobner basis can be constructed degree-wise. 

Definition 1. A finite set {gi, ...,gk} C R of homogeneous polynomials is a ho- 
mogeneous Grobner basis up to degree d of the ideal (gi,...,gk), if 

rem (S(gi, g/); gi, ...,g k ) = 

or deg (S(gi,gj)) > d, for all i,j = 1, 

Lemma 1. Let {gi, <j%} C R be a homogeneous Grobner basis up to degree d, and 
let p £ R be a homogeneous polynomial of degree at most d. Then, p £ (g%, ■■■,gk) 
if and only if Yem(p; g l7 ... 7 g k ) = 0. 



S{p,q) 



LCM(lm(p),lm{q)) 
W) 



P 



LCM(lm(p),lm(qj) 
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Proof. The paragraph preceding the definition implies that {gi, ...,gk} can be ex- 
tended to a Grobner basis Q of (gi,...,gk) by adding homogeneous polynomials 
whose degrees exceed d. Since deg(p) < d, we have rem(p; Q) = rem(p; g±, ...,gk). 
Since p S (Q) if and only if rem(p; Q) — by Theorem [2 the result follows. □ 

We see that in order to do Step (2) in the Basic Algorithm (or the correspond- 
ing steps in the Refined Algorithm) it suffices to know a homogeneous Grobner 
basis up to degree d of (P U Sd)- Our key theorem states that this Grobner basis 
can be constructed iteratively, as follows. 

Theorem 3. Let Q C R be a homogeneous Grobner basis up to degree d of (Q). Let 
p € R be a homogeneous polynomial of degree d, and p ^ (Q) . Then Q U {rem(p; Q)} 
is a homogeneous Grobner basis up to degree d of {Q U {p}}. 

Proof. Let r — rem(p;C?). Since p {Q) and all polynomials are homogeneous, we 
have r ^ 0, deg(r) = d, and (Q U {p}) = (Q U {r}>. 

By hypothesis, the S-polynomials of pairs of elements of Q are of degree > d or 
reduce to modulo Q. We now consider the S'-polynomials of r and elements of Q. 
Let g & Q. By definition of the remainder, we have lm[r) ^ lm(g). Therefore the 
^-polynomial of r and g is of degree > d — deg(r). Thus the claim follows. □ 

We obtain the 
New Algorithm 

(1) Compute m,d and a homogeneous Grobner basis Q of (P). Let Sd = LSd = 0. 

(2) For all power products b of degree d of elements of ISi U IS2 U • ■ • U LSd-i- 

(a) If rem(6; Q) <> then replace S d by S d U{b} and Q by 5U{rem(&; g)}. 

(b) If \Sd\ = rrid then break and return {Sd,ISd)- 

(3) For all b e B d : 

(a) If remO; Q) <> then replace S d by S d U {6}, by /Sd U {6} and 

gby£U{rem(&;£)}. 

(b) If ISdl = m d then break and return (S d ,IS d )- 

By Theorem [3] and induction, 5 is a homogeneous Grobner basis up to degree d of 
(P U £<$). Hence, in Step (2) (a) and (3) (a) one has rem(b;Q) <> if and only if 
b (P U 5<j). The New Algorithm is a dramatic improvement of the Refined 
Algorithm. However, in our examples this was still not enough. 

One should take more care in Step (2) of the New Algorithm. It simply says 
"For all power products b of degree d of elements of LS\ U IS2 U • • • U LSd-i" ■ Two 
questions arise: 

(1) How shall one generate the power products? 

(2) Is it necessary to generate all possible power products, or can one restrict 
the search? 

In very complex computations, the number of power products is gigantic. But 
usually only a small proportion of them will be eligible as secondary invariant. So, 
for saving computer's memory, it is advisable to generate the power products one 
after the other (or in small packages), rather than generating all power products at 
once; this answers Question (1). 

Apparently Question (2) was never addressed in the literature. However, it turns 
out that a careful choice of power products provides another dramatic improvement 
of the performance of the algorithm. Our choice is based on the following lemma. 
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This lemma seems to be well known, but to the best of the author's knowledge it 
did not appear in the literature and it was not used in implementations. 

Lemma 2. Assume that secondary invariants of degree < d are computed such 
that all of them are power products of irreducible secondary invariants. In the quest 
for reducible homogeneous secondary invariants of degree d, it suffices to consider 
power products of the form i ■ s, where i is a homogeneous irreducible secondary 
invariant of degree < d, and s is some secondary invariant of degree d — deg(z). 

Proof. Let p £ R be a power product of degree d of irreducible secondary invari- 
ants. Hence, it can be written asp = iq, with an irreducible homogeneous secondary 
invariant i of degree < d and some homogeneous G-invariant polynomial q of de- 
gree d — deg(?) (we do not use that q is a power product of irreducible secondary 
invariants). 

Recall that the secondary invariants generate the invariant ring as a {{P})— 
module. Hence one can rewrite q = qo + kiSi + ■ ■ - + ktSt, where qo € (P), k\, kt S 
K, and s\, ...St are homogeneous secondary invariants of degree deg(g). We obtain 
p = iqo + fci(zsi) + • • • + ktiist). Hence, rather than chosing p as a ((P))-modulc 
generator of R , we may choose is\, ...,ist, which, by induction, are all power 
products of irreducible secondary invariants. □ 

Improved New Algorithm 

(1) Compute nid- Let Q be a Grobner basis of (P). Let Sd = ISd = 0- 

(2) For all products b = i ■ s with i S I Si U • • • ISd-i and s 6 Sd-dcg(i) '■ 

(a) If rem(6; Q) ^ then replace Sd by Sd U {b} and Q by Q U {rem(6; Q)}. 

(b) If \Sd\ = md then break and return (Sd,ISd). 

(3) For all b e B d : 

(a) If remO; G) ^ then replace S d by S d U {6}, IS d by IS d U {b} and Q 
by £U{rem(&;£)}. 

(b) If \Sd\ — md then break and return (Sd,ISd)- 

This is the algorithm that is implemented as secondary_charO in the library 
f invar of Singular 3-0-2 [6], released in Juli 2006. In Step (2), the secondary 
invariant s may be a non-trivial powerproduct itself, hence, can be expressed as 
s = i s s', where i s is an irreducible secondary invariant and s' is (by induction) 
some other secondary invariant. Of course one should consider only one of the two 
products i s (is') and i(i s s') in the enumeration. 

Often one is only interested in the irreducible secondary invariants, which, to- 
gether with the primary invariants, generate the invariant ring as a sub-algebra. 
Therefore we implemented yet another version of the Improved New Algorithm 
in Singular 3-0-2, namely irred_secondary_charO. This algorithm computes 
irreducible secondary invariants, but does not explicitely compute the reducible 
secondary invariants. That works as follows. 

Let Qp be a Grobner basis of (P). In Step (2)(a) of the Improved New Algo- 
rithm, one replaces Sd by 5' c ;U{rem(&; Qp)}, rather than by S^Uffr}. In Step (3) (a) 
one replaces Sd by Sd U {rem(6; Qp)} and ISd by ISd U {b}. In the end, Sd does 
not contain secondary invariants, but normal forms of secondary invariants with 
respect to Qp. Since rem (rem(pi; Qp) ■ rem(p2j Qp)', Qp) — rem(pi ■ P2',Qp) and 
since a reduction modulo Q in Steps (2) (a) and (3) (a) also comprises a reduction 
modulo Qp, this maintains all informations that one needs for determining how 
many secondary invariants are reducible in Step (2) and for finding the irreducible 
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secondary invariants in Step (3). So in the end, ISd contains the irreducible sec- 
ondary invariants in degree d. This detail of our implementation very often saves 
much memory and computation time, as can be seen in Table [1] in Examples (1) 
and (6)~(9). In Example (8), we can compute the irreducible secondary invariants 
although the computation of all 31104 secondary invariants exceeds the resources. 

An example of Kemper (example (9) in the next Section) motivated us to further 
refine the implementation of the Improved New Algorithm. It concerns the 
generation of Bd'. If there are irreducible secondary invariants in rather high degrees 
d (in Kemper's example, there are two irreducible secondary invariants of degree 
9), it is advisable to generate not all of Bd at once, but in small portions. This will 
be part of release 3-0-3 of Singular. 

5. Benchmark Tests for the Computation of Invariant Rings 

5.1. The Test Examples. We already mentioned that some of our test examples 
arise in low-dimensional topology. This yields Examples (1), (7) and (8). For 
background information, see [12]. We will not go into details here, but just provide 
the matrices and primary invariants of our nine test examples. They are roughly 
ordered by increasing computation time. The ring variables are called Xi,X2, 
Let ei be the column vector with 1 in position i and otherwise. Our focus was 
not on the computation of primary invariants; note that in various examples the 
primary invariants are not optimal. 

(1) A 13-dimensional representation of the symmetric group S2 is given by the 
matrix 

M = (e 2 eiei3ei2eiie 8 eioe 6 e9e7e 5 e4e 3 ) 
Our primary invariants are 

Xg, £7 + £10, X 6 +X S , X 5 +Xn, X4 + X i2 , X3+X13, 
X1+X2, X 3 Xi 3 , X4X12, X S Xn, XjXiQ, X 6 X 8l X1X2 

There are 32 secondary invariants of maximal degree 6, among which are 
15 irreducible secondary invariants up to degree 2. 

(2) A 6-dimensional representation of S4 is given by the matrices 

Mi = (eie4e 5 e 2 e 3 e 6 ) 
M 2 = {e^eie^eaez) 

Our primary invariants are 

x-i +x 5 + x 6 , xi+x 2 +x 4l x 3 x 5 + x 3 x 6 + x 5 x 6 , 

X 3 X4 + X2X5 + X\X§, X1X2X4, x\ x\ + x\ x\ + X%x\ + X^X^Xq 

There are 12 secondary invariants of maximal degree 9, among which are 4 
irreducible secondary invariants of maximal degree 3. 

(3) A 6-dimcnsional representation of the alternating group A4 is given by the 
matrices 



Mi = 
M 2 = 



{e4e 1 e 5 e 2 e 6 e 3 ) 
{e2e 3 e\e§e4e$) 
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Our primary invariants are 

x\ + x 2 + x 3 + x 4 + x 5 + xq, x 3 x 4 + x 2 x 5 + x\Xq, 

X\X2 + X1X3 + X 2 X 3 + X1X4 + X2X 4 + X1X5 + X3X5 
+ X4X5 + X2X6 + X 3 X 6 + XiXQ + X5X6, 

X3X4 + x 3 xl + x\x b + x 2 x\ + x\x & + X\x\, 

X1X2X4 + X1X3X5 + X 2 X 3 X 6 + X4X 5 X 6 , 
X-yX<2 ~\~ X-^X^ ~\~ x^x^ ~\~ X-^X^ ~\~ x^x^ 

~\~ X^Xt^ ~\~ X^Xr^ ~\~ X-^Xt^ ~\~ X^Xq ~\~ X^Xq ~\~ X^Xq H - X^Xq 

There are 18 secondary invariants of maximal degree 11, among which are 
8 irreducible secondary invariants of maximal degree 5. 

(4) A 6-dimensional representation of the dihedral group D 6 is given by the 
matrices 

Mi = (e 6 e 5 e4e 3 e2ei) 
M 2 = (e 3 eie2e 6 e4e 5 ) 

Our primary invariants are the elementary symmetric polynomials. There 
are 120 secondary invariants of maximal degree 14, among which are 10 
irreducible secondary invariants of maximal degree 4. 

(5) A 8-dimensional representation of D s is given by the matrices 

Mi = {e 6 e 7 e & e b e 4 e 3 e 2 ei) 
M 2 = (e 4 eie2e 3 e$e 5 e 6 e 7 ) 

Our primary invariants are 

Xi + X 2 + X 3 + X4 + X5 + x e +x 7 + x», 

X4X5 + XlX 6 + x 2 x 7 + x 3 x s , x 3 x 5 + x 4 x e + XlX 7 + X 2 X S , 

x 2 x 5 + x 3 x 6 + x 4 x 7 + XiXs, xix 5 + x 2 x 6 + x 3 x 7 + x 4 x s , 

XlX 3 + X 2 X 4 + X 5 X 7 + X 6 X%, XlX 2 X 3 X 4 + x 5 x e x 7 x 8 , 
Xix\ + X2X\ + x\x 4 + X 3 x\ + x\x§ + XqX 7 + XyXg + X^Xg 

There are 64 secondary invariants of maximal degree 11, among which are 
24 irreducible secondary invariants of maximal degree 5. 

(6) A 7-dimensional representation of D\ 4 is given by the matrices 

Mi = {e 2 e 3 e 4 e 5 e 6 e 7 ei) 
M 2 = (eie 7 e 6 e 5 e 4 e 3 e 2 ) 

Our primary invariants are the elementary symmetric polynomials. There 
are 360 secondary invariants of maximal degree 18, among which are 19 
irreducible secondary invariants of maximal degree 7. 

(7) A 15-dimensional representation of S3 is given by the matrices 

Mi = (e2eie 3 e i e 7 ei4e 5 esenei 3 e d ei 5 e w e 6 ei 2 ) 
M 2 = (eie 3 e2e4e5e9e 8 e7e 6 ei 3 ei2eiieioei 5 ei4) 



10 



SIMON A. KING 



Our primary invariants are 

X1+X2+X3, XlX 2 + XlX 3 + x 2 x 3 , X1X2X3, 
2:10+2:13, 2:102:13, X 6 + Xg + Xu + X 12 + X 14 + X 15 , 
XllX 12 + X 6 X 14 + XgX 15 , XqXu + X 6 X 12 + 2:14X15, 
X 6 Xu + XgXi2 + XgXu + X12X14 + X 6 Xi 5 + X11X15, 



x^xgXi4 + X6X11X14 



XuX!2Xu 
J° 1 „,6 



X6X9X15 + X9X12X15 + 2:11X12X15, 
■x 7 ' 



X5X7 + X 5 X 8 + X 7 X 8 , X 5 X 7 X 8 

There are 1728 secondary invariants of maximal degree 17, among which 
are 76 irreducible secondary invariants of maximal degree 4. 

(8) A 18-dimcnsional representation of S3 is given by the matrices 

Mi = (e 2 eie3e4ei2eioe7eiiei4e 6 e 8 e 5 ei5e9ei3ei7ei 6 ei 8 ) 
M 2 = (eie3e2ei4e 8 e7e 6 e5e 9 eioei5ei3ei2e4eiiei 6 ei 8 ei 7 ) 
Our primary invariants are 
X1+X2+X3, xix 2 + X1X3 + x 2 x 3 , X1X2X3, 

X4+X9+X14, X4X9 + X4X14 + X9X14, X 4 XgX\ 4l 

X16+X17+X18, Xl 6 Xl7 + Xl 6 Xl 8 + X17X18, Xl 6 Xl 7 Xl 8 , 

X 6 +X 7 +Xl0, X 6 X 7 + X 6 Xl0 + X7X10, 

x 6 x 7 xi , 2:5 + X 8 + Xn + X12 + X13 + X15, 

X5X12 + X 8 Xl 3 + X11X15, X 8 Xn + X12X13 + X5X15, 

X5X11 + X 8 Xl2 + X5X13 + X11X13 + X 8 Xl 5 + X12X15, 

2:52:3X12 + X5X11X12 + X 5 X 8 Xl 3 + X11X12X15 + X 8 Xl 3 Xl 5 + X11X13X15, 

xt + xl + x e n +x e 12 +x e 13 + x 6 15 

There are 31104 secondary invariants of maximal degree 22, among which 
are 137 irreducible secondary invariants of maximal degree 4. 

(9) A 10-dimensional representation of S5 is given by the matrices 



Mi = 



/10 


{) 








0\ 





1 


1 

3 


1 

3 


1 

3 







1 

3 


2 
3 


2 

3 










2 
3 


1 

3 


2 
3 










2 
3 


2 
3 


1 

3 



















1 
















1 
















1 




\ 

\ 











1 / 













00100/ 




/I 












0\ 



1 

3 


2 
3 


2 

3 





2 
3 


1 

3 


2 
3 








2 

3 


2 
3 


1 

3 





1 


1 

3 


1 

3 


1 

3 

















-1-111 








{) 








-10 


1 





{) 








-1 10 



\ 











-10 

-1100/ 












M 2 = 



We are not listing the primary invariants here, as they are too big poly- 
nomials. There are 720 secondary invariants of maximal degree 22, among 
which are 46 irreducible secondary invariants of maximal degree 9. 
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Examples (2), (3) and (9) belong to a very interesting class of examples that 
was shown to us by G. Kemper [9]. For n € N, Let M n be the set of two-element 
subsets of {1, ...,n}. Then, one studies the obvious S n action on M n (or similarly, 
the obvious A n action), and one can try to compute the invariant ring Q[M n ) Sn 
(resp. Q[M„] A "). 

The 10-dimensional representation of S5 in Example (9) is a surprisingly chal- 
lenging problem. To simplify the computations, Kemper provided a decomposition 
of the representation into a direct sum of a 1-, a 4- and a 5-dimensional represen- 
tation. Without ad-hoc methods, the computation of secondary invariants for that 
problem has been beyond reach. The procedure (Irreducible) Secondarylnvari- 
ants of Magma V2.13-8 breaks immediately, since it requests 55.62 GB memory, 
while the memory limit of our computer is 16 GB. Our algorithm irred_seconda- 
ry_charO in Singular version 3-0-2 exceeds the limit of 16 GB while computing 
secondary invariants in degree 8. 

The total number of secondary invariants in Example (9) is not particularly 
large. The difficulties in Example (9) come from the fact that there are irreducible 
secondary invariants of rather high degrees. 

5.2. Comparison. We describe here how different algorithms perform on Exam- 
ples (1) up to (9). All computations had been done on a Linux x86_64 platform 
with two AMD Opteron 248 processors (2,2 GHz) and a memory limit of 16 GB. 
The computation of primary invariants is not part of our tests. Hence, in each 
example we use the same primary invariants for all considered implementations. 
We compare the following implementations: 

(1) secondary_charO as in Singular release 2-0-6. In Table [TJ we refer to it 
as "Singular (1998)". 

(2) secondary_charO as in Singular release 3-0-2, whith a small refinement. 
In Table [H we refer to it as "Singular (all sec.)". 

(3) irred_secondary_charO, as in Singular release 3-0-2, with a small re- 
finement. In Table [TJ we refer to it as "Singular (irr. sec.)". 

(4) Secondarylnvariants in MAGMA V2.13-8. 

Implementation (1) is due to A. Heydtmann [7] (1998) and has been part of 
Singular up to release 3-0-1. 

Implementations (2) and (3) are our implementations of the Improved New 
Algorithm explained in Section |4j They are part of Singular 3-0-2, released 
in Juli, 2006. Here, we test a slightly improved version, that saves memory when 
generating irreducible secondary invariants in high degrees. However, this only 
affects example (9); the performance in the other eight examples remains essentially 
the same, as the degrees of their irreducible secondary invariants are not high 
enough. 

Implementation (4) is due to A. Steel, based on [TD] or [H] or [3J. We consider here 
the MAGMA- version V2.13-8, released in October, 2006. There is also a function Ir- 
reducibleSecondarylnvariants in MAGMA, but computation time and memory 
consumption are essentially the same, in our examples. So, for the sake of simplicity, 
we do not provide separate timings for that function. 

Note that, after posting the first version of this manuscript, there was a new 
release of Magma containing an algorithm that G. Kemper developed in 2006. 
However, it seems that Kemper did not describe his algorithm in a paper yet. 
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Algorithm: 


C1 \ 


(21 


\°) 








Singular 


Singular 


Singular 


Magma 






(1998) 


fall ^pp i 






H/Api. 


(~\ \ 


0.55 s 


0.05 s 


0.03 s 


0.05 s 






8 62 MR 

U . VJ i_ 1\±. 1_J 


1.49 MR 


1.0 MR 


10 3 MR 


H/Xpl. 


\ Z ) 


0.05 s 


0.04 s 


0.04 s 


0.01 s 






99 MR 

VJ • \J \J 1 V 1 1 J 


96 MR 


97 MR 

\J m KJ 1 1 V 1 1 J 


7 05 MR 


H/Xpl. 


\ 6 ) 


0.48 s 


0.33 s 


0.3 s 


0.19 s 






2 97 MR 


1 95 MR 


1 96 MR 

_1_ . ij \J XV±J_J 


8 96 MR 


JzjXpi. 




6.55 s 


0.63 s 


0.32 s 


0.48 s 






12.29 MR 


2.47 MR 


2.97 MR 


9.09 MR 


Expl. 


(5) 


18.15 s 


10.53s 


9.69 s 


6.66 s 






45.79 MR 


10.61 MR 


17.0 MR 


31.82 MR 


Expl. 


(6) 


> 984 m 


100.4 s 


16.55 s 


118.51 s 






> 167 MR 


110.0 MR 


39.0 MR 


54.0 MR 


Expl. 


(7) 




268.9s 


20.94 s 


> 7 h 








872.7 MR 


35.1 MR 


> 15 GR 


Expl. 


(8) 




> 10 h 


50.7 m 










> 10 GR 


3.36 GR 


(259.5 GR) 


Expl. 


(9) 




6.42 h 


99.2 m 










10.74 GR 


7.35 GR 


(55.62 GR) 



Table 1 . Comparison of different implementations 



Meanwhile we implemented another, completely different algorithm in Singular. 
It will be part of Singular release 3-0-3 and often works much faster. E.g., it can 
compute Example (8) in 1.06 seconds. We describe this algorithm in [13] and also 
provide there comparative benchmarks using the new versions of Singular and 
Magma. 

Interestingly, in contrast to the corresponding Magma functions, irred_secon- 
dary_charO often works much faster and needs much less memory than seconda- 
ry_charO; see Examples (1) and (6)-(9). However, this is not always the case, as 
can be seen in Examples (4) and (5). 

In Table [I] " — " means that the computation fails since the process exceeds the 
memory limit; in examples (8) and (9), Magma requests the amount of memory 
that we indicate in round brackets. In some cases, we stopped the computation 
when it was clear that it takes too much time; this is indicated in the table by 
"> 

In conclusion, our benchmarks provide some evidence that the Improved New 
Algorithm has great advantages in the computation of invariant rings with many 
secondary invariants. Here, it marks a dramatic improvement compared with pre- 
vious algorithms in Singular or algorithms in Magma. In 3 of our 9 examples, 
it is the only algorithm that terminates in reasonable time with a memory limit 
of 16 GR. A particular benefit or our algorithm is that the computation of irre- 
ducible secondary invariants does not involve the explicit computation of reducible 
secondary invariants, which may save resources. 
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